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Abstract 

We give a complete classification of integrable potentials in the plane V = 
r^ 1 U{9) with U meromorphic and real analytic. In the more general case V 
only meromorphic, we make a classification of all integrable potentials pos- 
sessing a Darboux point c with V'(c) = — c and Sp(V 2 V(c)) C { — 1,0,2}, 
and then prove the non existence of other meromorphically integrable po- 
tential such that Sp(V 2 K(c)) c] — oo, 20]. An algorithm is given to improve 
this bound arbitrary, and allows to find algorithmically all integrable poten- 
tials in any set E of potentials posseding a generic property, the eigenvalue 
bounded property, which corresponds to the fact that for any potential in 
E, the lowest eigenvalue is always lower than some fixed explicit bound. We 
eventually present an analysis and some examples in the degenerate Darboux 
point case V\c) = 0. 

Keywords: Morals- Ramis theory, homogeneous potential, D-finiteness, 
higher order variational equations 

1. Introduction 

We want to find all meromorphic integrable homogeneous potentials of 
degree —1 in the plane. As we will not succeed to find all meromorphically 
integrable potential in this large class of potential, we will restrict ourselves to 
subsets of potentials satisfying a generic property, the "eigenvalue bounded" 
property. 
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The set of meromorphic potentials is the most general reasonable one, 
excluding potentials involving algebraic functions. It containts in particular 
all the rational potentials with among them the co linear 3 body problem, the 
co linear 4 body problem with symmetric masses (mi,m2, rri2, mi), etc... and 
all problems in higher dimension which possess some particular restriction in 
the phase space corresponding to a 2 degrees of freedom potential. We will 
then be able to solve all these problems at once, providing they satisfy the 
condition of "eigenvalue bounded" and that we compute this bound. For the 
following, we will only consider complex potentials and complex functions, 
except if specified otherwise. 

Definition 1. We will say that V is a homogeneous potential of degree —1, 
meromorphic in q\,q 2 ,r for r ^ 0, if 

V : g 2 , r) E C 2 x C* such that r 2 = q\ + g 2 } — > C 

and V is homogeneous of degree —1 in qi, q 2 , r. 

We use this definition because we want also consider the potential V 
which is invariant by rotation, and which is not a meromorphic potential. 
This is why we want to add the dependence in r, which will be needed by 
the way to consider polar coordinates. 

Theorem 1. Let V be a homogeneous potential of degree —1 meromorphic 
in q±,q 2 ,r for r / 0. Then it exists a 2n-periodic meromorphic function U 
such that for all (qi, q 2 ,r) G C 2 x C* ; r 2 = q\ + q\ we have 

V{.Q.i, <Z2j r ) = r U($) with rcos9 = qi, rsin9 = q 2 (1) 

Conversely, ifV can be written as ([T]), then V is a homogeneous potential of 
degree —1 meromorphic in qi,q 2 ,i r for r ^ 0. 

Remark that r, 9 correspond to the polar coordinates, and that for a given 
(<7i> <?2> r ) ; ® is on ly defined modulo 2n, which has no consequences because 
U is a 27r-periodic meromorphic function. 

Proof. Suppose that 

V(qi, q 2 , r) = r^ 1 U{9) with rcos9 — q%, rsin9 = q 2 , r 2 = q\ + q 2 
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with U meromorphic 27r-periodic. Then U can be written U(9) = f(expi9) 
with / meromorphic. So we replace then these by their expressions in Carte- 
sian coordinates. This gives 



V(qi,q 2 ,r)=r f 



qi + Hi 



As / is meromorphic, the function V is meromorphic in its 3 variables for 
any (qi, q 2 , r) G C 2 x C*. Conversely, we use homogeneity, and then we have 

V( qi ,q 2 ,r) =r-V(cos0,sin0,l) V(gi,g 2 ,r) G C 2 x C* 

using r cos 9 — qi, r sin 9 = q 2 . The function V(cos 9, sin 9, 1) is meromorphic 
27r-periodic. □ 

To a potential V we will associate a Hamiltonian system and the following 
differential equation system 

H(pi,p 2 ,qi,q2,r) = -(pl+pl) +V(qi,q 2 ,r) with r 2 = q\ + q\ 

q\ = Pi q 2 =Pi r = r~ l (piq x + p 2 q 2 ) 

9 TT -1 9 TT . 8 TT _] d TT 

Pi = ^T—H + r qi—H p x = —H + r q 2 ^-H 
oqi or oq 2 or 

The function H will be called the Hamiltonian, and this is a Hamiltonian 
system with 2 degrees of freedom. The function H is meromorphic for 
(Pi,P2,qi,q 2 ,r) G C 4 x C* and r 2 = q\ + q\. 

Thanks to Theorem [TJ it is sufficient and even more general to study the 
potentials of the form V = r^ 1 U{9) with U a 27r-periodic meromorphic func- 
tion. We can also ask why one should do this generalization. This is because 
the potential invariant by rotation is not a meromorphic potential because 
of the square root (in the distance to the origin). It is still meromorphic in 
qi, q 2 , r, and as we will see, it is very important that this potential belongs to 
the class of potentials we are investigating to do our "generic classification" . 
We call it only "generic" because it is useful only if some generic property 
of "eigenvalue bounding" is satisfied. This is linked to the fact that this 
classification rely on a conjecture (conjecture 2) that depend on some integer 
n and which can be proved algorithmically but only for a fixed n. 
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Definition 2. Let V be a homogeneous potential of degree —1 in the plane 
such that V = r^ 1 U(8) in polar coordinates with U 27r-periodic meromorphic 
function. We will say that V is meromorphically integrable if there exists a 
first integral I(qi,q2,Pi,P2,r) where / is meromorphic for (pi,P2,Qi,Q2,f) G 
C 4 x C* and r 2 = q\ + q\, and I is independent almost everywhere with H. 

In the following, we will often note ||<7|| 2 = qf + qi, even for complex values 
of q. In particular, this quantity is conserved by a rotation, even a complex 
one. 

Definition 3. We say that c = (ci, c 2 ) G C 2 is a Darboux point of V if 

d d 
— V(c) = aci — V[c) = ac 2 
dqi dq 2 

with a G C called multiplier. With homogeneity, we can always choose 
a — 0, —1. We say that c is non degenerated if a ^ 0. A homothetic orbit 
associated to c is given by 

qiit) = a<f>(t) Pi(t) = Ci<p{t) % = 1, 2 

with (j) satisfying the following differential equation 

The first order variational equation of H near a homothetic orbit is given by 

X{t) = ^VV(c)X(t) 

and after diagonalization (if possible) and variable change <p(t) — > t, the 
equation simplifies to 

2t 2 (Et + l)Xi - tXi = XtXi, Xi G Sp (V 2 V(c)) 

The main theorems of this article are the following: 

Theorem 2. Let V be a homogeneous potential of degree —1 in the plane 
such that V = r~ 1 U(9) with U a 2rr -periodic meromorphic function and 
C/(M) Cl IfV is meromorphically integrable, then 

v = - a em.* 

r 
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Theorem 3. Let V be a homogeneous potential of degree — 1 in the plane 
such that 

• V = -U(6) with U meromorphic 2n -periodic 

• 3c € C 2 such that cj + c\ ^ 0, V'(c) = -c and Sp(V 2 (V)(c)) C 
{-1,0,2} 

where Sp correspond the spectrum of a matrix. If V is meromorphically 
integrable, then V belongs after possibly rotation to one the following families 

V = - + - a , b 6 C, (a, b) ± (0, 0) V =- a e C* 
?i 92 r 

V = + + ? aeC*, e = ±l Eietannta 1987 0/ 

(gi + ezg 2 )' :i gi + e«g 2 

To prove such theorems, we will mainly use Morales Ramis theory about 
meromorphic integrability 

Theorem 4. (Morales, Ramis, Yoshida f^],f^],fj],f^j) If V is meromorphi- 
cally integrable, then the neutral component of the Galois group of the varia- 
tional equation near a non degenerated homothetic orbit associated to a Dar- 
boux point c with E — 1, c\ + c\ ^ over the base field C(0(t)) (where 4>{t) 
is the time parametrization of the homothetic orbit) is abelian at all orders. 
Moreover, if we fix the multiplier of the associated Darboux point to — 1, the 
Galois group of the first order variational equation has an abelian neutral 
component if and only if 

Sp (W(c)) C V^{k - 1)0 + 2), ke N j 

In fact, this is not exactly the same statement as the original theorem 
because we allow r to appear in the potential and in the first integrals, and 
because we precise explicitly the base field we use for Galois group compu- 
tations. There are two important arguments to use Morales Ramis proof 

• The first integral needs to have an expansion in series (or the quotient 
of two series in the meromorphic case) on the curve 

• The coefficients of this expansion will be functions of the time t, and 
the corresponding field will be the base field to be considered in Galois 
group computation. 
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The degenerate case will be studied in the last section. 

Proof. We know that c\ + c 2 . 7^ 0. The Morales Ramis theorem necessitate 
that the first integrals to be meromorphic along the curve. This is not the 
case but almost. In fact, a first integral I meromorphic in qi,q2,Pi,P2,r 
will be meromorphic along the curve except maybe at the point 0. The 
parametrization of the curve is 

qi(t) = Ci<f)(t) pi(t) = Ci<p(t) i = l,2 



We remark in particular, that on this curve, we have r = \Jc\ + c|0(£) (this 
square root correspond to a choice on which branch of the potential the 
homothetic orbit is, because the potential is bi-valuated as a "function" of 
qi, 52)- Moreover, on the surface r 2 = q\ + q\ we have 

d qi r = ^ d q2 r = ^ (2) 

Here we see that the formulas ([2]) are always well defined on the curve be- 
cause a/ c\ + c\ 7^ 0. So, the initial form of / along this curve will still have 
coefficients meromorphic in (0, <fi) (with possibly an essential singularity at 
= 0,oo), because such an initial form is computed by doing derivations 
and then evaluation on the curve. We now need to pass from meromorphic 
coefficients to rational coefficients. This is because all computations we are 
able to do are over a rational base field, and not a meromorphic one. 

Let us note J the initial form of /. After the variable change <f>(t) — > t, 
the normal variational equation of order 1 becomes of the form 

2t\l + t)y" -ty'-Xy (3) 

where A is a parameter corresponding to an eigenvalue of the Hessian matrix 
V 2 V(c). The function J is a homogeneous function in y, y, and its coefficients 
are in meromorphic in t, a/1 + t" 1 for t 7^ 0, 00. We remark that this square 
root does not appear in the coefficients of the variational equation. So just 
by multiplicating J by itself changing the valuation of the square root, we 
can suppose that J is meromorphic in t for t 7^ 0, 00. All singularities of 
equation ()3]) are regular. So the Galois group over C(t) of equation ([3]) is 
exactly the Zariski closure of the monodromy group. If we remove small disks 
D , .Dqo on the singularities, we will have that the first integral J will 
be meromorphic on 

C \ (D U U Doo) (4) 
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If we restrict ourselves to take paths in this set instead of taking path in C 
to compute the monodromy group, this will not change anything (because 
we have chosen the disks small enough, and then all curves in C \ {0, —1, 00} 
are homotopic to curves in the set (jlj)). If the Galois group is SX 2 (C), then 
the monodromy group is non commutative, and then such a first integral J 
cannot exists. 

In their theorem, Morales, Ramis use the Kimura table [gj , which in fact 
correspond to Galois group computation over the base field C(t) of equa- 
tion ([3]). They found that the Galois group of equation Q is SX 2 (C) except 
for the specific values of A in Theorem H] (for which the Galois group is 
abelian). 

□ 

We will find in particular that a non degenerated Darboux point of V 
corresponds to some #0 G [0, 2tt[ such that 

U'(6 Q ) = and U(0 ) ^ 

One of the main difficulties comes from the fact that function U can be arbi- 
trary flat near 9q, and we then absolutely need higher variational equations at 
an arbitrary order. Indeed, previous classification of integrable homogeneous 
potentials in the plane for the degrees 3, 4 have been done for polynomial po- 
tentials 0,0, and this corresponds to a finite number of parameters. Here 
the function U can be anything, and thus we need an infinity of integrability 
constraints to prove these theorems. 

Example 

n2 n 

This potential is integrable at order 2 near all Darboux points (using in- 
tegrability table in |9j). Moreover, near 9 = (which corresponds to a 
Darboux point), this potential is integrable at least at order 2n — 1. The 
analysis of the other Darboux point (9 = 71) is even more difficult because 
Sp(V 2 (V)(c)) C {2, (2k -l)(k + 1)}, integrable at order 2, and integrability 
at order 3 has an explicit dependence on k. We will prove that it is in fact 
non integrable at order An — 3 near 9 = 0. For Theorem |2l we do not have as 
hypothesis the existence of a Darboux point. Because of this, we will need 
the following theorem 
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Theorem 5. Let V be a meromorphic homogeneous potential in the plane of 
degree — 1 such that V = r~ 1 U(9) with U meromorphic 2rc -periodic and real 
analytic. Then it exists 9q such that 

u(6 )^o u\e ) = o ^^<o 

Proof. We have that U(R) C R and U is The function U is periodic, 
then it exists a minimum and a maximum for U . Suppose first that U is not 
constant. Then maxC/ > mint/. We have 3 cases 

• maxC/ > mm.U > 0. Then we pose 9 such that U(9 ) = maxU 

• maxU > > minU. If maxU ^ 0, we pose 9 such that U(6 ) = 
max £7, otherwise we pose #o such that U(9 ) = minU 

• > max U > min U. We pose then 9 such that U(9 ) = min U 

Knowing that maxC/ > minf/, we get U(9 ) ^ 0. Then in all cases, we have 

U"%) 
U(9 ) 

Knowing that is an extremum, we get 

U(9 )^0 U'(9 ) = ^y<0 

which gives the theorem. □ 

Remark 1. We have then that the first theorem is inside the second one. 
In particular, the computation in polar coordinates gives 

If V is meromorphically integrable, then the eigenvalues at Darboux points 
should belong to 

| (p-i)o»+') , peN | = { _ 1|0[2[5|9|14 ,... } 



< 
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Knowing that 

U"(6 ) 
— — < 

u(o ) - 

we get the condition using Morales Ramis integrability condition H] 

Sp(VV(c)) = {2,-l} 
This case is in the hypotheses of the second theorem. 

2. Eigenvalue Bounding 
We will denote 

M = {V(r,6) = r- 1 U(9) with U meromorphic 2n periodic} 

Let V G M.. We note d(V) the set of Darboux points c of V with multi- 
plier —1 and such that ||c|| 2 ^ 0. Let c G d(V), then we have Sp (V 2 V(c)) = 
{2, A} and we note 

n _ f A if A G R 
1 — oo otherwise 

We consider E C At a subset of .M and we define the following 
A(E) = sup inf A(c) 

We say that the problem of finding all meromorphically integrable potentials 
in E is a bounded eigenvalue problem if 

A(E) < oo 

We have A(A^) = oo because of the following family 

V(r, 6) = r" 1 ((1 + a) - 2ae w + ae 2id ) , aGl, 

for which there exists only one Darboux point c = (1,0) and the correspond- 
ing eigenvalue is A = 2a — 1. 

Remark 2. Theorem [5] has then proved the following 

A ({V = r- x U{6) U ^ meromorphic 2vr-periodic U(R) C R}) = -1 



Here we want to present the following way to prove non integrability. 
First, you find a majoration of A(E) of your problem. Then you run a 
program to check non existence of integrable potentials with eigenvalues up 
to this majoration (this program and associated conjecture are in the section 
"The other eigenvalues"). Finally you look if inside the family E, there is a 
potential given by Theorem [2J You do not even have to check Morales Ramis 
condition of order 1 given by Theorem |H Thanks to this a lot of integrability 
problems can be proved. Moreover, this eigenvalue bounding property is 
more or less generic. This is because the Maciejewski-Pryzbylska relation on 
eigenvalues of Darboux points holds under generic assumption 1(| (even for 



non polynomial potentials), and this relation implies eigenvalue bounding. 
This is why we can call Theorem [3] and the partial proof of conjecture 2 
used together as a "generic" classification, because it holds for all eigenvalue 
bounded potentials sets. 



3. Non degeneracy of higher variational equations 

We will first recall some properties of the solutions of the first order 
variational equations. After diagonalization and in the integrable case, the 
equation is the following (after fixing the energy E — 1) 

2t 2 {l + t)y" -ty' - -(n - l){n + 2)y neN (5) 
We make then the change of variables 

v / TTF T — >t 

and the differential equation becomes 

- l)y + 2ty - l -(n - l)(n + 2)y = 

A basis of solutions is given by (P n , Q n ) where P n are polynomials in t and 
the functions Q n are 

QM = pjt) J ^ _ ^ (ty it 

The functions Q n are multivalued except for n = which will be a special 
case. Indeed, the Galois group of (j5J) in this case is Id instead of C. 
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The polynomials P n can be computed using the "Rodrigues" type formula 

1 d n - 1 



which gives a normalization for the leading term of P n that we will choose 
now and the functions Q n can be written as 

Q n (t) = e n P n (t) arctanh (\) + Vl "' ' ' 



J J t 2 -l 

with W n being polynomials, and e n a real sequence given by 

4T n n(n + 1) 



-n — 19 



Lemma 6. (already proved in fM]) Let F e C(zi) [22] and 

f(t) = F \ t, arctanh 
We consider the following field extension and monodromy group 

K = C (t, arctanh 0J , J fdtj , G = a(K, C(t)) 
J/G zs abelian, then 

Res F ( t, arctanh ( — ) + a ) — Va G 



da t=co \ \t 
where Res correspond to the residue. 

Proof. We will consider two paths, the "eight" path <j\ around the singu- 
larities —1 and 1, and the path 02 around infinity. At infinity, the function 
F (t, arctanh Qj + a) will have a series expansion of the kind 

J F ^t, arctanh ^jj + c^j dt = a n (a)t n + r(a) In t 



because the function arctanh ( |J has a regular point at infinity. Let us now 
consider the monodromy commutator 
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We have that a?* (f) = F (t, arctanh (}) + 0) and a 2 (\nt) = hit + 2in. We 
deduce that 

<r(f) = f + rtf)-r(P) 

This r(a) corresponds to the residue of F (t, arctanh (|) + aj at infinity. If 
the monodromy is commutative, then the commutator a should act trivially 
on /. This is the case only if r(j3) — r(0) V/3 G Z. The function r is a 
polynomial in (3, so r(j3) — r(0), V/3 G C. This gives us the formula of the 
lemma 



□ 



—Res F t, arctanh - + a = VaGC 

CjQ; t=oo V V f 



Definition 4. Let V = r 1 U(9) be a meromorphic homogeneous potential 
of degree —1 such that c = (1, 0) is a Darboux point with multiplier —1, and 
we note Sp(V 2 \/(c)) = {2, l/2(i - l)(i + 2)}. Let k G N* be fixed and (VE k ) 
the k-th order variational equation near c. We suppose (VE^-i) integrable. 
Then we will say that the integrability constraint of (VEk) is non degenerated 
if 

^-Res (t 2 - l)\Q % + e % aPi) k+1 ^ 
oa t=oo 



Now we can ask why such definitions. This is because non degeneracy im- 
plies a rigidity property. If we take two integrable potentials "close" enough 
(meaning that enough derivative on some Darboux point are equals), then 
they should be equal. Remark that our non degeneracy constraint could be 
generalized in other cases. 

Behind this, there is the following philosophy. Let V : C n — > C be 
a homogeneous potential such that c = (1, 0, . . . , 0) be a non degenerated 
Darboux point. Let k G N* be fixed and (VE^) the fc-th order variational 
equation near c. We suppose (VE^-i) integrable. Then we will say that the 
integrability constraint of (VEk) is non degenerated if the abelianity of the 
identity component of the Galois group of (VE k ) depend explicitly on the 
k + 1-th derivatives of V on c. This notion of "explit dependence" depend 
on the degree and dimension you consider, and also the rigidity result you 
want. Such a procedure in fact suppose that we want to use integrability 
conditions to express the k + 1-th derivatives with lower order derivatives, 
and then apply a recurrence property. Once first derivatives are determined, 
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the higher ones will be determined every time this non degeneracy condition 
is satisfied (and in practice, it is generically satisfied if the first order Galois 
group is not too simple). 

Theorem 7. Let V : C 2 — > C be a integrable meromorphic homogeneous 
potential of degree — 1 such that c = (1, 0) be a Darboux point with multiplier 
— 1. Suppose it exists k such that integrability constraint of iVE^) is non 
degenerated Wk > k Q . Let V be another integrable meromorphic homogeneous 
potential of degree — 1 such that 



Qi+j 

dq\dq 3 2 



V{c) 



Qi+j 

dq\dq 3 2 



V(c) Vi + j<k 



Then V = V 

Proof. Let us prove it by recurrence. Suppose that all (VEj) are integrable 
up to i = k — 1 > ko — 1 and non degeneracy of the integrability condition 
at all orders > ko. The k-th variational equation is written (one of them) 



X(t) 



2 
A 



X(t) 



( i JS) vi yi-j \ 



2^ i|/|.!2 
i=2 Vy j=0 



\ 



E 



1 



i=2 Vy j=0 



i-3 



J 



(6) 



Here we use the following abusive notation. The functions X^X 1 ^ 3 are not 
unknowns of the equation. They corresponds to solutions of linear differ- 
entials equations systems, with possibly non homogeneous terms. A more 
complete description of the construction of higher variational equation sys- 



tems can be found in [ll|]. We also suppose by recurrence that the coefficients 



,0') 



< k — 1 are already determined. These coefficients corresponds to 
the derivatives c^c^V^c), i + j < k. Let us now prove that the coefficients 



a 



b): are uniquely determined. 
We have that the constants a 



(j) lO') 



k j u k+i are on ^y in equation (Q. They do 
not appear in the linear differentials equations systems whose solutions are 
the X^X 1 ^ 3 . This is because they correspond to the term of highest order in 
the series expansion. Then the monodromy of (VE^) is an affine function in 
a k ' i^k ' We now use Euler relation 



qid qi V + q 2 d q2 V = -V 
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and we derive it in q% or q 2 enough, and then evaluate it on (g 1; q 2 ) = (1, 0). 
This gives us all derivatives of order < k + 1 except 

ak+i 

Then all coefficients are known except b^ . We need to determine it using 
the fact that V is integrable. We write the solution of the second equation 
of © 

X(t) = X hom (t) + X partl (t) + 6 fe ° 'X part2 (t) 

and we apply to it a monodromy commutator a = o~f "cr^Vi* ^ (with the 
same notation as in the proof of Lemma [6]). This gives 

o{X) = a(X hom ) + a(X partl ) + bf } a(X part2 ) 

Now let us look at X part2 . It can be computed and one solution is 



X part2 (t) = J(t 2 -l) k Q k+l dt 



We now apply lemma [6] which says that the monodromy element a add to 
such function the constant 

G(a) = Res (t 2 - l) k (Qi + e ia Pi) k+1 - Res (t 2 - l) k Q h+1 

t=oo t=oo 

This is not zero by hypothesis and so we get this equation 

= <y(X) — X = a(X hom ) — X hom + a(X partl ) — X partl + bf^G(a) 

that we can solve in bf\ So b® is uniquely determined. The meromor- 
phic property of V allows us to conclude. All derivatives on c are uniquely 
determined, then V is uniquely determined. □ 

This rigidity allows to produce uniqueness theorems, but never existence 
theorems because we never know if all integrability conditions are compati- 
ble or not. For classification results as Theorem El there is then a "praying" 
phase, where it is needed to find explicit integrable potentials in the litera- 
ture, whose first derivatives can be arbitrary. Direct search help a lot like 
[H, but if it is not successful, all previous work is almost useless (because 
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only unicity is proved, not existence). This suggests to carefully read the 
literature before beginning to prove similar theorems as |3] by using the non 
degeneracy condition [7J This phase cannot be algorithmically made because 
computing higher variational equation will only give a series expansion of the 
corresponding potential (and because computer limitations, with often less 
than 10 terms). 

3.1. Application to the eigenvalue 

Theorem 8. Let V be a meromorphic homogeneous potential of degree —1 
in the plane. Suppose that c = (1,0) is a Darboux point with multiplier 
— 1, Sp(V 2 V(c)) = {2,0} and that V is meromorphically integrable. Then 
V = l/ gi . 

Proof. We just need to prove the non degeneracy theorem. The functions 
Pi, Qi for the eigenvalue are the following 

Pi = 1 Qi — arctanh 

We need to look at the following residue 

Res {t 2 - l) k+1 (arctanh f - ] + a - 

t=oo y yt ) t 2 — l 

and it should be independent of a. The easiest coefficient to study (and non 
trivial) seems to be the coefficient in a k+l . We find after simplification 

S k = Res (t 2 - l) h+1 f(k + 2) arctanh 

By expanding, we remark that the second term always give a zero residue. 
Indeed, in the expansion, the faction simplifies and we get a polynomial. 
Then, we will compute 

i 

S k = Res (k + 2) (t 2 - l) k+1 arctanh (-) = — f {f _ dt > 
t=oo \t J 2 J 

-l 

The last equality is produced with the expansion of arctanh Q) at infinity. 
We deduce that 

S k ^ VA: > 1 



G- 



kt + It 
t 2 - 1 
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So we now know that there is a unique potential with (1, 0) as Darboux 
point with multiplier —1 and eigenvalue 0. The potential l/q\ satisfy these 
conditions, and is integrable because it is invariant by translation. □ 

Conclusion. So, after rotation, an integrable potential V with a zero 
eigenvalue near a non degenerate Darboux point correspond to the potential 

V=-, aeC* 
Qi 

3.2. Application to the eigenvalue 2 

In the case of the eigenvalue 2, the Hessian matrix should be diagonaliz- 
able. This is a condition for integrability of first order variational equation. 
First we will prove the non degeneracy condition. 

Theorem 9. For eigenvalue 2, the integrability condition is non degenerate 
at order k > 3, and degenerated at order k = 2. 

Proof. We need to look at the following residue 

Res (t 2 - l) k+1 [ + 6ta + 6t arctanh - 

t=oo \ t 2 -1 \t 

and this residue should be independent of a to prove non degenerescence. 
We will look at 

= Res (t 2 - l) k +H k+1 (-^fr 4 + 6t arctanh 

which corresponds to the coefficient in a k+1 (after simplifying a non zero 
factor). But this sequence is not always non zero. We will also need to look 
at another one 

= Res {t 2 - l) k+ H k (-%^ + 6* arctanh (-X) 

t=oo \ t 2 — I \t J J 

Then, we want to prove 

S« ^ or S { p ^0 \/k>2 
More precisely, we will prove that 

S<£ ± and SgVi ^0 VA; > 1 

16 



We use Mgfun to build a recurrence and then an explicit expression of these 
sequences. We find the following recurrence for 

64(2n + 3)(2n + l)(6n + ll)(n + 1) 2 /W~ 
(20736n 5 + 152064n 4 + 439200n 3 + 622752n 2 - 431784n - 116328) 
/(ra + 1) + 36(6n + 5)(3n + 5) (3n + 4) (6n + 13) (6n + 17) /(n + 2) 

This recurrence can be solved explicitly and gives the formula 

(1) vrr(2n + 2)27-" 

^2n — 



72T (n+l)T(n+f) 



This expression never vanish. We do the same for S^ +1 . We find a 3-th 
order recurrence and solve it 

( 2) Tr27' n T (2n + 3) ^ ( (3k + 4) r (k + 5/3) T (k + 7/3) \ 

2 " +1 ~ 3456r(n+|)r(n+f) ^ \ (* + 1) (* + 2) (2k + 3)T (k + f)T (k + ^) ) 

Using this expression, we find that S^n+i never vanish for n > 1. This prove 
the non degeneracy condition for order > 3. At order 2, the two formulas 
vanish. Since in this case the residue is a polynomial in a of degree at most 
2, this implies that the residue is constant. So the a derivative is zero and 
integrability condition is degenerated. □ 

Now we know that after fixing the 3-th order derivatives of the potential 
V, the integrability condition completely determine the potential. Still we 
need to find such potentials. We already know the potential 

1 , o,06L 

which is integrable. Computation gives that Darboux points have the eigen- 
value 2. So we need to prove that after rotation, all possible 3-th order 
derivatives can be produced. As we will see, this is only "almost" true. There 
is an unexpected set of parameters impossible to produce. These parameters 
correspond in fact to the Hietarinta potential. 

Theorem 10. Let V be a meromorphic homogeneous potential of degree —1 
in the plane with c = (1,0) a Darboux point of V with multiplier —1 and 
V 2 y(c) = 27 2 - Then it corresponds to a potential of the form 

- + -, a,beC* (7) 

Ql 12 
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after rotation except if V possess the following expansion 

V(c + q) = V(c) + agi + (qf + q\) + dq\ + 3dq iq j ± 2irfg| + o(q 3 ) 
Proof. We expand V r onc=(l,0) which gives 

V{c + q) = V(c) - qi + (qf + q\) + aq 3 + bq\q 2 + cq.qj + dq 3 + o(q 3 ) 
We use Euler equation 

qid qi V + q 2 d q2 V = -V 
and after derivation, we have 

2d hl V + qi d lthl V + q 2 d 1 , 1 , 2 V = -d ljX V 
2d 1<2 V + qidi tlt2 V + q 2 d lt2<2 V = -d 1>2 V 
2d 2)2 V + q\di^ 2 V + ^2^2,2,2^ = -d 2 , 2 V 
Thanks to this, we have the following values 

dx,i,iV(c) = -6 d 1X2 V(c) = d h2>2 V(c) = -6 
Then the series expansion of V on c is always of the form 

V(c + q) = 1 - gi + {q\ + q\) - {q\ + 3 gi g 2 2 + dql) + o(q 3 ) 

where d G C We want now prove that such an expansion can correspond to 
the expansion of the potential ((?]) after rotation. So we will make a rotation 
of the coordinates qi,q 2 - After rotation, the potentials (j7|) can be written 

+ , w 2 + z 2 = l,a,&GC* 

The condition of possessing a Darboux point at c = (1, 0) with multiplier — 1 
implies that this family of potentials can be written 

z 3 W 3 99 

+ T-5 ^ r i W 2 + Z 2 = 1 



(z 2 + w 2 )(zs — wy) (z 2 + w 2 )(wx + zy) 
We make series expansion near c = (1,0) and by identification, we get 

-Z 2 +W 2 2 2 

= a w + z = 1 

This produce the solution 




4 + d 2 + V4rf 2 + d 4 



A + d 2 

which is valid for d 7^ ±2z. For d = ±2z, there are no solutions. □ 
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3.3. The irreducible case 

There is still one open possibility. It corresponds in fact to the Hietarinta 
potential 

V = ^M±^ + ^_ aeC*,e = ±l 

[qi + eiq 2 ) 6 qi + eiq 2 

Contrary to the other ones, this one can not be decoupled after rotation. 
It still possess an additional first integral quadratic in momenta, but the 
quadratic term is (pi + ip2) 2 , and effectively, this quadratic form can never 
be decoupled after rotation (any rotation, even a complex one). This is why 
it is a special case. 

Theorem 11. Let V be a meromorphic homogeneous potential of degree — 1 
in the plane with the following expansion at c = (1, 0) 

V(c + q) = V(c) + b qi + {ql + q\) + dq\ + Mq x ql ± 2idq\ + o(q 3 ) 

If V is meromorphically integrable, then it corresponds after rotation to a 
potential of the form 

V = + ql l + — aGC*, e = ±l 

{qi + cm) Qi + mi 

Proof. We already know that if such potential satisfy these conditions, it is 
unique (after fixing the sign of e). Let us check that the Hietarinta potential 
satisfied them. 

V = a- ; — r^ + a — e = ±1 

[qi + 6iq 2 y qi + eiq 2 

We will study the case e = +1. After rotation, we get 

V = a-, — + ao , o,oet 

(qi+m) qi+iq2 

The condition of having a Darboux point at c = (1,0) with multiplier — 1 
gives 

y_ 1 l! + <ll , * 



2 (qi + iq 2 ) 3 2(qi + iq 2 ) 

We compute the series expansion at c = (1,0) and this gives exactly the good 
expansion. □ 
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4. Case of the eigenvalue —1 

The case of eigenvalue —1 is much more difficult, because higher vari- 
ational equations do not seem to have particular properties, and the non 
degenerescence theorem does not apply. In fact, we will even prove the re- 
verse, that the integrability condition is always degenerated. We need to use 
a completely different method. We already guess that this case will corre- 
spond to the potential invariant by rotation. Then, to see something special, 
it seems to be a good idea to compute all these higher variational equations 
in polar coordinates, because it corresponds to the action-angles coordinates 
for this potential. This will be a priori much more difficult, but in fact, the 
eigenvalue —1 will give us lots of simplifications. 

Theorem 12. Let V = r~ 1 U(8) be a meromorphic potential homogeneous of 
degree —I in the plane with c = (1,0) a Darboux point of V with multiplier 
— 1. Suppose that Sp(V 2 V(c)) = {2,-1}. IfV is integrable, then 

V = r~ l 

Proof. The potential V = r~ 1 U(9) gives the following differential equations 
in polar coordinates 

f - r 2 = ~U(6), rO + 2f9 = ^U'(9). 

Let us linearize this equation near a homothetic orbit corresponding to the 
critical point of U. We suppose moreover that U"(0) = 0, which corre- 
sponds to Sp(V 2 V(c)) = {2, -1}. We get at first order 

We choose the normalization of c such that the multiplier is —1 and we make 
the variable change (pit) — > t which gives 

2 + l)f-ir = ^r, 2(t + l)0+ + 4^ = 0. 

Then we make the change of variables y/1 + t _1 — > t which gives 
\ (t 2 - l) f + 2ft - 2 r = 0, \ (t 2 -1)9 = 0. 
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Of course these equations are integrable (because they correspond to an 
integrable case of Theorem H] and the solutions are 

r(t) = C x P 2 + C 2 Q 2 , 6{t) = C 3 t + C 4 

Now we take a look at second-order variational equations. We find 



- (t 2 - l) f + 2tr - ^9 2 = 2 r - 3(t 2 - l)r 2 

i§ + I(t 2 - lye + 2tr6 + (t 2 - iye = ^-^-jU {3 \o)e 2 



(8) 



by using the notation abuse consisting to consider that the terms of degree 2 
in r and 9 are replaced by the solutions of the first-order variational equations. 
The first equation integrates because 

2 2 3 

corresponds to add to the solution the constant — |C| and that the mon- 
odromy of 

J {t 2 - l)Q\dt 
is commutative. For the second equation, we find 

^9+2t(C 1 P 2 +C 2 Q 2 )C 3 +(t 2 -l)(C 1 P 2 +C 2 Q 2 )C 3 = ^-^U^(0)(C 3 t+C 4 ) 2 
The solution can be written as 

2 J j -2t(C 1 P 2 +C' 2 Q 2 )C 3 -(^-l)(C 1 P 2 +C 2 Q 2 )g3+^ (3) (0) ( °%* 
We have that 

1 3+2 



P 2 = 4t Q 2 = -t arctanh - + % 8 

The terms in P 2 are polynomials and integrate well. For the terms in Q 2 , we 
find the following expression 

J J 2tQ 2 + {t 2 - l)Q 2 dtdt = J J jj(3t 2 - 1) arctanh ( jj - ^tdtdt = 

— (t 2 - l) 2 arctanh ( - ) - — t 3 
32 v ; \t J 32 
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Then the only term left is 

For we integrate and we find that 

J I l^(3)(0) (C3 ^+^ 4)2 rftdt ECt, arctanh f , In (t 2 - l) (9) 

Then the second-order variational equation is always integrable, and more- 
over we have that 

• If U^(0) ^ then the Galois group of © is C 2 

• If U^{0) = then the Galois group of © is C 

We have in fact already more or less proved the following 

Theorem 13. Let V be a meromorphic potential homogeneous of degree —1 
in the plane with c = (1,0) a Darboux point ofV with multiplier — 1. Suppose 
that Sp(V 2 1/(c)) = {2, -1}. Suppose that U®(0) = 0Vi = l.Jfe toen toe fact 
that the identity component of the Galois group is abelian or not does not 
depend on the value ofU^ k+l \Q). 

Remark 3. This imply the degenerescence of the integrability condition at 
all orders. 

Proof. The case of order 2 correspond to the previous proof. Let us look 
now at order k. We pick in the equations the non homogeneous terms where 
U^ k+1 '(0) appear. The only equation where such term appear is the following 

l-._l ^(0) gfc 
2 2 t 2 - 1 

(we have suppressed all other non homogeneous terms for which U^ k+1 \0) 
do not appear). We conclude with the following property 

J J dtdt G c t*> arctanh ^jj , In (t 2 - l)] 

which can be checked using recursive integration by parts. This term has then 
a commutative monodromy, and then the fact that the identity component of 
the Galois group is abelian or not does not depend on the value of U^ k+1 \0). 

□ 
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Remark 4. Still we remark that in fact this integral will grow a little the 
Picard Vessiot extension, and the Galois group become at least C 2 for a 
generic value of [A fe+1 )(0), but this does not give us the non integrability. 
This is precisely why this case is particularly difficult. We cannot use non 
degenerescence properties, and so we need to keep these unknown coefficients 
and go higher in the order of variational equations. 

Let us look now at order 3. We find 
and yi satisfies the following system of differential equations 

yi = 2y 2 V2 = + 77^77 2/3 = 2 v ^ _^ 

with z = £/ (3) (0) and by using the notation abuse consisting to consider that 
the terms of degree 3 in r and 9 are replaced by the solutions of the first 
order variational equations. We will only consider a "normal perturbation" , 
which consists to put r = in all terms of degree 3. Moreover, to simplify, 
we will choose 9(t) = 1 as solution for the first order variational equation. 
We also have that the terms 

^(t 2 - l)r9 + 2tr9 + (t 2 - l)r9 

will have no impact on the Galois group because these terms are solutions of 
a linear differential system of equations with non-homogeneous terms, but 
these non- homogeneous terms vanish for r = 0, 9 = 0. Then the solutions of 
these systems will be solutions of the second order variational equation, and 
we have already determined that the corresponding Galois group is abelian. 
Because of that, we come down to the study of the equation 

The term in 9 3 integrates like ([9]) and so this term will not give a non-abelian 
Galois group. The terms left are 

1-- 1 z z9 3 z9 2 9 
7;° = ^ — 7?/i y\ = 2y 2 2/2 = 2/3 + 75 m = 2- 



2 2t 2 -T * * * *° t 2 -l *° t 2 
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We solve and we find the following solution (a well chosen one) 

yi (t)=2zjj2j Ht) t2 ^ 6 i {t) dt + ^dtdt 9{t) = JI ¥ ^y 1 (t)dtdt 

The 9{t) in terms of degree 3 is replaced by the solution of the first order 
variational equation, and we had chosen 9{t) = 1. We have then 

yx(t) = 2z ^-tarctanh - i In (t 2 - l)J 

We know that f ^-ryi{t)dt is in the Picard-Vessiot field and we have 

J — Uarctanh + l - In (t 2 - l) J dt = 

^ln 2 (t- 1) -ln21n(t-l) - dilog(t + 1)- 
Iln 2 (t+l)-lln(t-l)ln(t + l) 

All the terms are in C[t, arctanh (|) , In [t 2 — 1)] except one, the dilogarithmic 
term 

dilog (*+!)= [ Ht + 1) dt 



t 

With the same idea as in Lemma [61 we get that this term has a noncommu- 
tative monodromy because of the following residue in 

^ ln(t + 1) + a 
ties = a 

t=o t 

which depends explicitly in a. 

Here we have then proved that the 3-th order integrability condition is 
[/( 3 )(0) = 0. We remark this is a very special case because the 4-th order 
derivatives do not appear in the constraint. We will now prove that a very 
similar scheme appear at higher order. 

Let us prove now Theorem 1121 by recurrence. Suppose that £/«(()) = 
Vz = l..k and let us look at the 2k — 1-th order variational equation. We 
want to prove that f/^ +1 ^(0) = 0. We have already initialized this recurrence 
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with k = 2. In this equation, there are the derivatives U^(0) for i — k + 
1 . . . 2k. We have already proved that the coefficient f/( 2fc )(0) will have no 
impact on integrability. We get then (ignoring terms with no influence on 
integrability, and considering only the normal perturbation) 

We use the abusing notation as before and so the 9 k+l are solutions of differ- 
ential equations system with non homogeneous terms coming from lower or- 
ders. As before, we will take 9(t) = 1 as solution for the first order variational 
equation, and thus we will have 9 = 0. This produce a huge simplification. 
We can choose for example the solution for the terms 

-{t 2 - l)r9 + 2tr9 + {f - l)r9 

which are solutions of homogeneous linear differential systems. We then only 
have the following 

1- f l(^(0) t+1 
2 d -2^k\ t*-l 9 

Still, these 9 k+l are not easy to compute a priori. We build the differential 
system for these 9 k+l by posing 

yf = QOQk+i-j j = 0...k + i 

and we have then 

k o jj(k+l+n) (r\\ 

yf = (k + i — j)yf +l + j ^ i L^ ^f 1 

n=0 

The degree in 9 of the second term is 2k + n + % — 1. But we study the 
variational equation at order 2k — 1, then we can suppress all terms of order 
higher than 2 k — 1. Then there is only one term in the sum 

yf = (k-j)yn i+ jl^^9^9^ yf = (k-3)y% i > 1 
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The only non homogeneous terms are of maximal degree 2k — 1, and we had 
chosen 6(t) = 1 as solution for first order variational equation. Then all term 
with 9 vanish. We have then only 

y ^ = + ______ y)' ^{k-jjy)^ (i,j)^(0,l) 

The system of equation number i, i ^ for the yj l \ is a linear homogeneous 
differential equation system. We can then take the following solution 

y f = Vi ^ 1 

We look now for % — 0. We have 

?/o 0) = M 0) 

yf } = (k - j)yfl j — 2. . .k 
We find the following solution of this system (very well chosen) 

,(°) 



2 f/( fc + 1 )(0) 



and we remplace it in the initial equation 



M(t arct anh + I In (t 2 - l)) 



4f/ (fc+i)( ) i / A\ I \ 



After solving it, we get that 

[/(fe+i)(0) 



^/^(tarctanhQ+iln^-l)) 



(fc-1)! 

should belong to the Picard-Vessiot field. This is the same integral as be- 
fore, and its monodromy is not commutative. Then a necessary integrability 
condition is that £/( fe+1 )(0) = 0, which completes the recurrence. □ 

The function U is meromorphic, we have that U^ k \0) = VA; G N* and so 
this implies that U is constant. Then 

V = r- 1 
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5. The other eigenvalues 

To find all integrable potentials, we would need to study all the other pos- 
sible eigenvalues. But in the other cases, there are no special simplifications, 
but moreover, we do not know integrable potentials with these eigenvalues. 
Then we can suppose that such potentials do not exists, and we can also 
make a stronger supposition that at some order, there is a condition impossi- 
ble to satisfy. First of all, we will give some application theorems for the non 
degenerescence property. This property cannot be checked for all possibilities 
for the moment, but with an axi-symmetric assumption, such property can 
be proved. 

Theorem 14. Let V = -U(0) be a homogeneous potential of degree —1 in 
the plane in polar coordinates with U meromorphic 2n -periodic. We suppose 
that 

• V is symmetric with respect to the straight line generated by 9 = 

• U(0) ^ 

// V is meromorphically integrable, then V belong to a discrete family 
v - a ° v - ai 1 v - ° 2 cos ^ V V V 

vq - — ,vi — 7^7; V2 77^7, V 3 , Va, v 5 , ... 

r r cos(0) r cos(2#) 

More precisely, V = aVi, a e C* with 




Remark 5. We say nothing about existence of the potentials Vi (for % > 3). 
We only prove that for each eigenvalue \{k — l)(k + 2), k e N, it exists at 
most one axi-symmetric integrable potential. 

Proof. We just need to prove the non degenerescence property for odd orders. 
Indeed, for even orders, we use Euler relation which gives us all derivatives 
except one, the derivative in the normal direction to the straight line 6 = 0. 
But for the variational equation of even order k, this maximal order derivative 
is then of odd order k + 1. This derivative is then automatically because 
we suppose the symmetry. The non degenerescence is written 

|-Res (t 2 -l) k (Q n + ae n P n ) k+1 ^0 
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We look at coefficient a k of the residue. 



e k n (k + l)Res (t 2 - l) k Q n P k = \e k n {k + 1) / P k+1 dt 

t=oo Z J 

-1 

by doing the Taylor expansion of arctanh (|) at infinity and recognizing that 
this sum can be written as this integral. The integer k is odd, the polynomials 
P n are never identically 0, then this coefficient never vanish. This prove non 
degenerescence. □ 

Conjecture 1. Let V = -U(9) be a homogeneous potential of degree —1 in 
the plane in polar coordinates with U meromorphic lix-periodic. We suppose 
it exists 9o such that U'{9q) = and 

U "( 9 °) le(-(n-l)(n + 2),^N,n>3 



U(9 ) " 1 2 

Then it exists a k e N* such that the variational equation at order k is not 
integrable. 

Conjecture 2. Let V = -U{9) be a homogeneous potential of degree —1 in 
the plane in polar coordinates with U meromorphic 2tt -periodic. We suppose 
it exists 9o such that U'(9 ) = and If 

l ^M-l= 1 -(n-l)(n + 2) noddn>3 

then the 5-th variational equation is not integrable. If 

U"(9 ) 1 

— 1 = -(n — l)(n + 2) n even n > 4 
U(9 ) 2 

then the 7-th variational equation is not integrable. 

The second conjecture seems to be extremely difficult. Solving it would 
classify completely integrable homogeneous potentials of degree —1 in the 
plane, and would probably allow with some generalization for other degrees 
to close completely the search of integrable homogeneous potentials (with at 
least some assumption on Darboux points). A partial proof up to the 5-th 
variational equation would lead to classification of axi-symmetric integrable 
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potentials. This would also lead to numerous theorems in higher dimen- 
sion for potentials posseding discrete symmetry groups. We will prove these 
conjectures for n = 3,4, 5, 6, and give all corresponding expansions of V in- 
tegral) le at order k (with k < 4 or < 6 depending on parity). This proves 
by the way the nonexistence of the potentials V3, V4, V5, V§ of Theorem [TU 
Moreover, we present an algorithm to prove this conjecture for an arbitrary 
fixed number of eigenvalues. 

Proof. (Partial proof of conjectures 1 and 2) We will suppose that c = (1, 0) 
is a Darboux point and that V has been normalized such that the multiplier 
of c is —1. We note Sp(V 2 ^(c)) = {2, A}. The corresponding angle in polar 
coordinates is 6 = 0. We begin by odd index eigenvalues, because we will 
just have to look at order 5. This correspond to A = 5, 14. We know that the 
3-th derivative in q<i should be zero for integrability at order 2. The Euler 
relation for V gives the following expansion near c 



where z is an unknown coefficient. This is this coefficient that our program 
will determine. To make more readable the results, we will write the solution 
with an expansion in polar coordinates (even if the variational equations are 
computed in Cartesian coordinates). We have for the moment (let us recall 
that the k-th order variational equation gives constraints on the k + 1-th 
derivatives) 



We want now to compute variational equations to very high order to see 
if we can always be integrable at high order by choosing good values for 
the derivatives. Let us build a general algorithm for computation of higher 
variational equation. We already know that we will need to compute the 
differential systems associated to 



This produce linear differential equations in yij t k,i with i + j + k + I < n 



V=- qi + -{2ql + \ql) - q \ - -\ qi q 2 2 + qf + 3\q 2 iq 2 2 + z q \ + o(g 4 ) 




Ui,j,k,l — Xi X2 X 1 X 2 
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where n is the order of the expansion. We use the following formula 

, , 4:ityi,j,k,i 4jtyij,k,i, 

Ui,j,k,l — KVi+l,j,k-l,l -r iyi t j + i,k,l-l f2_i 1 2 — 1 

i add(add(coeff(coeff(W^i,g™),gf)^i iiifc+ni i +p ,n = O..oo),p = O..oo) + 
j add(add(coeff(coeff(W^ 2 ,g^),gf)|/ iii _ lifc+nii+ p,n = O..oo),p = O..00) 

where 

W 1 = subs (q! = (t 2 - l)q u q 2 = (t 2 - l)q 2 , ^-^ d qi V 

W 2 = subs = (t 2 - l)q u q 2 = (t 2 - l)q 2 , ^-^ d q2 V^J 

To simplify, we look at the system with % + k and j + I fixed. The Vij,k,i i n 
the equations have alH + k and j + I equal or higher. This produce a block 
triangularization of our system. 

Practical computation 

We already know the solutions for i + j + k + I = n. We can then take 

Vi,j,k,i = QpQ l p if i = k = 

Vi,j,k,i = otherwise 

This choice correspond to a normal perturbation, with as much as possible 
functions Q (this is because they are more complicated, so we can expect 
more integrability conditions from them). By linearity, we also can make the 
even better choice 

Ui,j,k,i = coeff((Q p + e p aP p y (Q p + e p aP p ) 1 , a m ) if i = k = 

Ui,j,k,i = otherwise 

which select in advance the coefficient in a in the final residue to compute. 
This allows to simplify the computation by doing it separately. Another 
difficulty is that even if the system should always have a solvable Galois 
group, it is very difficult to solve. We will then just search solutions of the 
form 

T 



Vi,j,k,i e C(f) 



arctanh 
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with the degree in arctanh less than m, except for the equation in i = k = 
and j + l = 1. This is in fact a very reasonable hypothesis. Behind this, there 
is the supposition that the Galois group will not grow before becoming non 
abelian. This hypothesis is reasonable, because if the Galois group grows, 
then some logarithm should appear at infinity. We have already seen that this 
phenomenon is exceptional in [9|. We also suppose that the degree in arctanh 
does not grow. The question of integrability at order n arrives only when 
solving the last equation. This last equation is solved with standart method 
of variation of the constant, and the integrability condition is established 
using Lemma |6j Let us precise that all these conditions are sufficient to 
assure the termination of the algorithm, but are in no way necessary. In 
all cases, if the algorithm gives a solution, the solution is always valid. One 
important precision is that we need to use only one solution (although well 
chosen) . 

We find the following expansions 

and that V3, V5 are never integrable at order 5. To prove this, the easiest way 
is to consider the case m = 1, and m = 3, and this give two incompatible 
constraints for the 6-th derivative. 

In the case of the even index, we have the following expansions at first 
order 

V A = -(1 + 56 2 + b6 3 + o(9 3 )) v 6 = -(l + —6 2 + b9 3 + o(9 3 )) 
r \ 1 r \ 2 J 

We have here a free parameter b. This is the main reason of the greater 
difficulty When we arrive at order 5, we also get two incompatible conditions 
for the 6-th derivative, except for some well chosen b. We find then the 
following expansions (for r = 1) 

1 oW 363467 , 112035, 

— d&Vi = b 6 H b 

5! 9 4824000 8576 
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with 



l 9 ^ 



216926052083 l4 279352141289, 2 
-b A + b 2 



10224685080000 



54531653760 



4715685295 
24563808 



R4b) = 



158469311 
97702546320000 



372429603 
868467078400 



45927 
2729312 



= 



68250852673 ,, 98831601, 
:6 + b 



4257725150000 



3475694 



6! 
with 



10915637473609903 
5190230823727250000 



6605928379884787, 2 19638863047783 
-b 2 + 



1271076936423000 



10039110960 



_ 198715111646995383 4 

6 ^ ^ ~~ 2772435940648535262500000 + 
1448561702310687 2 270431334600 

7921245544710100750 + 3128145145507 ~ 

So for each expansions for eigenvalues 9, 20 there are exactly 4 possible 
expansion integrable at order 5. We remark that this set of possible expan- 
sions is invariant by the transformation 6 — > —9 which is expected because 
this transformation is a symmetry and does not change integrability or non 
integrability of the potential. For higher orders, we need then to suppose 
that Ri(b) = 0. The explicit expression of b is too complicated to do effi- 
ciently the computations at higher order. So, we first remark that b appear 
only in the non homogeneous part of the linear differentials systems we have 
to solve. This means that all our computations can be done in 



Q(6) 



1 , f 1 

t, — , arctanh - 

' t 2 - 1 ' \t 



X,t, 



— .arctanh (i 

- 1 \t 



im*)) 



if Ri is irreducible (which in our case has always been true but it can of course 
be tested for each eigenvalue). We then just have to write for each differential 
system the non homogeneous part as a linear combination of 1, b, b 2 , b 3 , then 
solve it for each term and then sum. Indeed, the solutions of each term will 
be in (at least we assume it with our assumption on Galois groups) 



1 , ( 1 

t, -= , arctanh - 

' t 2 - 1 It 
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When we reinject the solution in the next differential system, there could be 
multiplications by b, so we reduce the non homogeneous part modulo Ri. We 
find at the end 

1 7 _ 57826741017348283 25932696791821703 &3 
7! 9 4 ~ 893377392990720 + 100504956711456000 



1 d 7y _ H8828154548524498748866853503827777 & 

7! 9 6 ~ 431797756299715943933989778480280 + 
8633140425176867273801758981735627411 ,„ 



52895225146715203131913747863834300000 



There is only one condition at order 6 and it can be satisfied. At order 7, we 
have again two constraints, and this time we do not have a free parameter to 
fix anymore. These two conditions cannot be simultaneously satisfied. □ 

Classification of Galois groups 

Order: 1 2 3 4 5 6 7 



A odd ^ c ^ 



A even ^ —1, 2 £■ 



A = 



A = 2 



C 



,n. a. n.a. n.a. 

C -C *C^n.a. 

^^n.a^^n.a. 
r C * C * C — - n.a. 

n.a. n.a. n.a. n.a. 

k C C C C * C — ► n.a. 

w n.a. n.a. n.a. n.a. n.a. n.a. 
- C - C C C C 

n.a. n.a. n.a. n.a. n.a. 

-c — *c — -c — -c — -c — *c 



33 




The graph is built using the following rules: At each order k, if there is 
branching, the branch choice depend on some condition on the k + 1-th 
derivatives or lower. The branch noted n.a. correspond to non abelian Ga- 
lois groups. The tree is infinite for eigenvalues —1, 0, 2, but at the end, only 
the lower branch correspond to integrable potentials at all orders (and these 
are truly integrable). 

Remark. - These expansions of the potentials are unique et allow for each 
given potential to precisely compute the order at which it is integrable (here 
in the case of the eigenvalues 5,9,14,20). This classification is then only 
proved for the eigenvalues —1, 0, 2, 5, 9, 14, 20, and for bigger eigenvalues, this 
corresponds the conjecture 2. But the algorithm we used can compute the 
same for an arbitrary fixed eigenvalue. The proof of conjecture 2 is then pos- 
sible for many eigenvalues, even if the computation cost is great. This leads 
to a practical way to study in all generality eigenvalue bounded problems, at 
least in the case where the bounding is effective (known explicitly). 

6. Other type of potentials 

We have here only written about non degenerated Darboux points. Now 
we will take a look at the degenerate case. 

Theorem 15. Let V = \U{6) be a meromorphic homogeneous potential of 
degree —1 in the plane in polar coordinates with U meromorphic 2n -periodic. 
We suppose it exists 9q such that 

U(6 ) = U'(6 ) = 

If V possess a first integral I, meromorphic for (pi,P2, Qi, ?2, r) G C 5 , r 2 = 
Qi + <?2 an d independent almost everywhere with H , then U"{0q) = 0. Con- 
versely, if U"(6q) = 0, then the identity component of the Galois group of 
variational equation near 9 is abelian at any order. 
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The fact that the first integral should be meromorphic for r = is in fact 
a serious limitation. It is possible that the potential itself is not meromorphic 
for r = 0, for example when U is not a rational function in exp i9. So in this 
sense we cannot call this a "reasonnable" definition of integrability, because 
the class of possible new first integral does not contain the Hamiltonian. 
So it seems that such a definition is reasonnable only when U is a rational 
function in expz#, which corresponds to V rational in g 1; g 2 ,r. Remark that 
the proof cannot be so easy as Theorem H] because the variational equation is 
not Fuchsian, and then we cannot use the trick with the monodromy group. 

Proof. Suppose first that U"(9q) ^ 0. The first order variational equation is 
the following 

Xi = X 2 = — jj^-X 2 

The first equation is clearly integrable. For the second one, we make a linear 
variable change and this gives 

v=$v (io) 

This equation is not Fuchsian because is not a regular singularity. Using 
the Kovacic algorithm is in this case a priori not enough. Indeed, the Ga- 
lois group over the field C(t) will detect only rational first integral. Let us 
prove now that if (ITUj) admit a first integral of the form I(y,y) with / an 
homogeneous rational fraction with meromorphic coefficients in t, then / has 
rational coefficients in t. 

Infinity is a regular singularity of ( flQj) and the Frobenius expansion near 
infinity gives (after variable change t — > ~) 



v (t) = c 2 ln t ( - + — t + —t 2 ) +c 2 ( - - — t 1 —t 2 ] + 

y U 2K ' V 4 32 768 y V* 64 2304 J 

ci ( - + — t + —t 2 ) + o(t 2 ) 
\4 32 768 J 

There is a logarithmic singularity. Suppose it exists a first integral meromor- 
phic in t and homogeneous in y,y of (flQl) . After variable change t — > |, 
there could be an essential singularity in (but not multivaluation) in the 
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coefficients of I. We pose 



J = t 2 ( - + —t + -^—t 2 + o(t 2 )) y - t 2 ( 1 + -t + — t 2 + o(t 2 ) ) y = 
V 8 96 3072 / V § 192 / 

The function J is the only first integral of degree 1 whose coefficients are 
univaluated in 0. Indeed, if there was another one (independent), then y, y 
would be univaluated for all choices of C\,C2, which is not true because In t 
is infinitely multivaluated. Let us look now at higher degrees. We consider 
I(y,y) a first integral and yo a solution such that I(yo,yo) = 0. Then 2/0/2/0 
is an algebraic combination of the coefficients, so 2/0/2/0 is at most finitely 
multivaluated, so this implies that 

Vo = Vafi a^o 

Then we can factorize in / 

i(y, y) = (t 2 yi,oy - t 2 yi, y)I(y, y) 

This strictly decrease the degree of the numerator. We can also do the same 
procedure with the denominator by considering I~ l instead of /. We conclude 
by recurrence that 

I = AI$ A G C*, keZ 

Then the coefficients of / are meromorphic at infinity (so they are rational). 
The Galois group of f lTU]) is SX 2 (C), and then it do not exists such first 
integrals (non constant). Then V is not meromorphically integrable. 

Now suppose the reverse, that U"(6q) = 0. The first order variational 
equation is written 

x x = x 2 = 

We now use the lemma 

Lemma 16. The algebra C[t, |, In t] is stable by integration. 

Proof. We consider / G C[t, |,ln t] and we write it as a linear combination 
of terms of the type 

rin(t) m n G Z, m G N 
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If n > 0, then we use integration by parts to decrease m until 0. If n < 0, 
We use integration by part to increase n up to n — — 1. We then have the 
formula 

Then all functions in C[t, |, In t] have a primitive in C[t, |, In t]. □ 

To conclude, we remark the following phenomenon: The solutions of 
higher variational equation are in fact solutions of non homogeneous lin- 
ear differential equations and the non homogeneous terms are produced only 
using products of lower order solutions and functions t~ k . So the solutions 
always live in some algebra in which we take recursively integrations. So 
we apply the method of variation of constants to find the solutions. More- 
over, the wronskian of X 2 = equal to 1 (and also for the higher variational 
equations matrices), then we never have to divide. So all solutions belong to 
C[t, |,ln t] which is stable by integration. Then the Picard Vessiot field is 

Ki = C(t) or C(t, In t) d = id or C 

The Galois group is virtually abelian at all order. □ 

Remark 6. This type of proof appears to be very general. Indeed, if some 
potential appears to be integrable at all order near a particular solution, 
without known first integral, the Picard Vessiot field is often not growing. 
If the coefficients of the potential are not well adjusted to avoid creating 
further monodromy, it is probably because it is not possible. The solutions of 
higher variational equations are in fact solutions of non homogeneous linear 
differential equations and the non homogeneous terms are produced only 
using products. So the solutions always live in some algebra in which we take 
recursively integrations. For homogeneous potentials and non degenerate 
Darboux point, the algebra is given by the following process 



Ao = C 



t. 



t 2 - 1 



Ai+i = I Aidt 



where f Aidt D Ai is the algebra generated by all primitivations of functions 
in A { . We have then in particular that the Picard Vessiot of (VEi) (in the 
case of all the eigenvalues belong to Morales Ramis table) is always contained 
in the fraction field of A^. These algebras containts in particular all the 
polylogarithms functions that give integrability constraints (however, it is 
not clear that Ai is only generated by polylogarithms). 
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7. Conclusion 

Eventually using this analysis, all eigenvalue bounded problems can be 
algorithmically studied. The two conjectures represents the last open ques- 
tions about homogeneous potentials in the plane of degree — 1. They could 
in principle be proved using the D-finiteness property of the functions P n , Q n 
(the fact that they satisfy linear polynomial recurrence and differential equa- 
tions), but in practice direct computation seems to be way out of reach for 
the moment. No counter example have been found to conjectures 1, 2, and 
numerical computations strongly suggest that the integrability conditions at 
order 5 and 7 are never compatible. With these theorems we see that very 
often, the difficulty is not the number of parameters, but specific family of 
potentials (with even as low as 1 parameters) which exhibit the non eigen- 
value bounding phenomenon. 

Some examples of potentials integrable at all order near all Dar- 
boux points 

In the case of non degenerated Darboux points, if we admit conjecture 2, 
it will not be possible to find non integrable potentials which are integrable 
all orders. We then need to find homogeneous potentials either having no 
Darboux points at all, either only multiple degenerated Darboux points (the 
second derivative should vanish). The functions U(0) — F (e* e ) with 

F{z) = h(z n ) h fractional linear, n e N* 

/oz l 
-dz < i < j - 2, n e N* a e C* 
yz oty 

have no critical points. The functions U{6) = F (e* e ) with 

F{z) = h((z n - a) m ) - h(0) h fractional linear, m > 3, n G N* a G C* 

have only degenerated critical points satisfying the integrability constraint. 

These examples show that there are still open questions about integrabil- 
ity, but the difficulties do not rely on Morales Ramis theory but on the search 
of Darboux points. Potentials without non degenerate Darboux points are 
not common, but they still exist and a complete classification of them seems 
to be difficult. 
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